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Abstract 

Quantum Monte Carlo simulations are used to investigate the two- 
dimensional superfluid properties of the hard-core boson model, which show 
a strong dependence on particle density and disorder. We obtain further 
evidence that a half-filled clean system becomes superfluid via a finite tem- 
perature Kosterlitz-Thouless transition. The relationship between low tem- 
perature superfluid density and particle density is symmetric and appears 
parabolic about the half filling point. Disorder appears to break the super- 
fluid phase up into two distinct localized states, depending on the particle 
density. We find that these results strongly correlate with the results of sev- 
eral experiments on high- Tc superconductors. 
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Quantum systems of interacting bosons in clean or disordered media have been used to 
study a variety of interesting problems, including the superfluidity of ^He, superconductor- 
insulator transitions in thin films, and vortex dynamics in type-II superconductors. Recently, 
the relation between superconducting transition temperature (Tc) and carrier density [6) 
in various copper-oxide superconductors was found to obey a universal bell-shape curve 
This finding also shows that the superfiuid-like condensation of charged bosons (lo- 
cal Cooper-pairs) is crucial for high-Tc superconductivity, based on some common features 
of copper oxides: the extremely short coherence length implying point-like Cooper pairs, 
and the layered structure which confines the carriers mainly to two-dimensional (2D) Cu02 
layers. From the analysis of /iSR and transport measurements, Schneider and Keller p| 
argued that such condensation of extreme type-II superconductors belongs to the classi- 
cal XY universality class in three dimensions, and Zhang and Sato surmised that the 
bosons involved are actually bipolarons. In this model. Cooper-pair bosons whose binding 
energy is much larger than condensation energy are assumed to exist below Tc(onset), and 
superconductivity appears via the phase coherence of pre-existing bosons. 

Here, within the boson framework, we propose a quantum XY universality class in two 
dimensions. We examine a 2D lattice boson Hubbard model specified by hard-core repul- 
sive interactions and random potentials. The phase diagram depends on the temperature, 
particle density, and also on the degree of disorder. First and foremost, when the particle 
density is varied, the phase diagram should be perfectly symmetric about the half filling 
level, because of the particle-hole symmetry. Second, the hard-core bose gas has mathe- 
matical analogy to the quantum spin problem 0]. From previous studies of spin-1/2 XY 
model one may expect that, for a clean boson system, the Kosterlitz-Thouless (KT) 



transition [|T0| takes place and 2D superfluidity arises below a finite critical temperature 
Tkt- Also, regarding different bosonic systems, there are early numerical works suggesting 
a KT transition in a Coulomb gas model and in a soft-core model [|T^. Third, at suf- 
ficiently strong disorder, dense bose gas is expected to exhibit a localized phase called the 
Bose glass Through a quantum Monte Carlo (MC) simulation, this paper provides 



further support for the superfluid KT transition of the hard-core boson model in the half fill- 
ing case, and sheds light on density-modulation and localization effects in 2D superfluidity. 
Our simulation results strongly correlate with recent experimental data in copper-oxide su- 
perconductors, including the Tc vs. 6 relationship and the nature of superconducting phase 
transition. 

The Hamiltonian of hard-core bosons is expressed as 

t ^ 

-H"boson = -- + + '^i^i + -^HC (l) 

for a square lattice of size LxL. Here the first sum is the kinetic energy of bosons, where 
the hopping constant t = h'^/2mla'^ with effective mass and lattice spacing a. The 
second part represents the potential energy from onsite disorders with a uniform distribution 
V G [—A, A]. The hard-core interaction i^HC inhibits the double occupancy of bosons at 
each site. This interacting bose gas can be transformed into an equivalent spin-1/2 XY 
system in a randomly varying magnetic field, such that i^boson -f^spin = ^tJ2{i.j){S^Sj + 
SfS]) + Ei ViSi, by setting = + iSy, a = - iS'^ , and n = + 1/2 with the spin-1/2 
operator S = (5^ 5^ 5") §. 

To simulate the equilibrium state of such a qunatum system efficiently, we use the path- 
integral approach based on the Suzuki- Trotter transformation The 2D bosonic sys- 
tem (under periodic boundary conditions) is transformed to a 3D (torus with space and 
imaginary-time dimensions) classical system of boson world lines. Possible paths of the 
world lines are topologically characterized by particle number and winding number, where 
the winding number {W) is defined by counting how many times the world lines wind around 
the torus in space directions. For the microcanonical ensemble of bosons (fixed particle den- 
sity rife), we carried out MC sampling of the world lines, by utilizing an algorithm previously 
developed for the spin model. ( See Ref. [0 for detail of the simulation method, and Ref. ||T^ 



for some technical remarks, respectively.) The superfluid density Us corresponds to the free 
energy change due to twisting the phase of order parameter along one lattice boundary, and 
this quantity can be computed directly from winding number fluctuations 



n, = {k^T/t){W^) (2) 



in the path-integral representation |15|. For the MC samphng, together with conventional 



local moves of world lines, we also made global update of winding number W |11,12], such 
that all the transitions to different values of W are allowed. In contrast, most previous spin 
simulations [^,|^ were limited to a fixed condition of = while the particle number was 
variable. Since the present VT-variable algorithm is ergodic regarding winding number, we 
can evaluate rig exactly for a finite-size lattice. This computational approach is complemen- 
tary to earlier quantum MC studies and exact diagonalization study [|l3l on the similar 
model. 

Figure 1 shows the numerical results of a clean half-filled system (A = 0, rib = 1/2) with 
different size lattices (L = 4,6,8). The temperature dependence of 2D superfiuid density 
is shown in Fig. 1(a). For T < Tkt = 0.40(5) t/ZcB? ns{T) grows remarkably and becomes 
almost independent of the lattice size, while, for T > Tkt, it is suppressed with increasing 
size. At low temperatures, ng reaches 0.27(1) for the L = 8 lattice, compared to 0.28 given by 
the exact diagonalization of the much smaller {L = 4) lattice |]13[. Thus, even at the ground 



state, the superfiuid fraction is at most ng/rih ~ 0.54 in the thermodynamic limit, and the 
normal-fluid component remains large. This indicates a strong renormalization effect due 
to quantum fluctuations, while the thermal fluctuations cause the vortex-antivortex pairs to 
unbind at high temperatures. The value of Tkt is quite consistent with the KT universal 
jump condition, which satisfies ^^(Tkt) = '^k^Ti^iT / nt denoted as a dotted line in Fig. 1(a). 
Moreover, to get a clearer sign of the KT transition, we calculated the temperature derivative 
of Us using the fluctuation formula 

di/3ns)/d/3 = {W^){E)-{W^E) (3) 

where E is the total energy and (3 = l/k^T. As seen in Fig. 1(c), the peak value of d{[3ns) / d[3 
grows rapidly near T = Tkt, as the lattice size increases. In sufficiently large lattices, Eq. (3) 
reduces to the energy difference /S.E between periodic and antiperiodic boundary conditions. 
It should be noted that, for the classical XY model, an analogously divergent sign of /S.E 



was used as a decisive evidence for the KT transition |T^. On the other hand, the peak 
value of specific heat tends to saturate at C^, = 0.62(1) just above Tkt with the increase in 
size, as shown in Fig. 1(c). Although it is well known that free bose gas in two dimensions 
never condenses at finite temperatures, these numerical observations definitely support the 
view that the hard-core nature of interactions between bosons leads to the superfluid phase 
through the finite temperature KT transition, while the quantum fluctuations are rather 
significant at low temperatures. 

Figure 2 shows how the superfluid density depends on the bose-particle density and 
the degree of disorder (A), at low temperature T = 0.25t//cB- Note that the relationship 
between the superfluid density rig and the doping density n?, is approximated well, although 
not exactly, by the parabolic form oc (rij, — no)(l — Uf, — Uq) with some constant Hq. 
Here Uq = 0.06, 0.13, 0.13, and 0.15 for various degree of disorder A = 0, 1.0, 1.25, and 
1.75, respectively. In comparison, in a one dimensional (ID) pure system, the exact result 



Ug = sm(7rniy) / TT is obtained [|T^. The symmetry about half filling is rigorous in both clean 
and disordered systems, since the model involves the particle-hole symmetry due to hard- 
core interactions. In the presence of disorder the parabola seems to be preserved, and in 
addition the onset of superfluidity can be found at a common value of uq ^ 0.1, despite 
large reduction of maximum Ug at the half filling. Since the transition temperature Tkt 
is generally proportional to ns{T = 0) in the 2D system, the above parabolic property 
implies that the relation between and particle density takes a universal form, which is 
accompanied with a common offset density in the disordered case. 

The disorder always localizes bosons to suppress the superfluidity. Nevertheless the 
disorder-induced localization appears to differ qualitatively between dense and dilute bose 
systems in two dimensions. Figure 3 shows the dependence of on disorder (A) at a fixed 
low temperature. The superfluidity of the moderately dense (half-filled) system is insensitive 
to weak disorder. At stronger disorder A > 2.5, the superfluidity disappears, but the size 
dependence is very weak even when the disorder is rather strong. This result is consitent 



with earlier studies [[I3i|l8[ which support the existense of the Bose galss phase for dense 



systems. In contrast, the dilute system at 1/8-filling {rih = 1/8) is rather sensitive even 
to weak disorder. Here is rapidly suppressed as the lattice size is increased. This size 
dependence implies that the correlation length of superfluidity is smaller than the lattice 
size L = 8. The insets of the flgure show the spatial distribution of boson density in 
each case. For the dilute region in the vicinity of Uf, = no, most bosons appear to form a 
superfluid cluster, as depicted in the upper inset. The size of the cluster is comparable to the 
localization length C,i of the Anderson transition. For instance, using C,i = a/ log(A/t) in the 
limit of strong disorder, C,i ~ 5a for A = 1.25. This observation indicates that the 2D weak 
localization is dominant at densities below uq. On the other hand, for a half-flUed system 
with strong disorder, such clusters assemble into a percolative network, as shown in the lower 
inset. This behavior is supposed to be associated with critical fluctuations in the vicinity 
of the Bose-glass transition point. Here, when the glass correlation length, larger than the 
Anderson's localization-length, exceeds the lattice size, the superfluidity disappears. By 
varying the particle density under the hard-core condition, we have thus obtained a possible 
tuning from Anderson glass to Bose glass in two dimensions. It is interesting that the 
similar tuning was previously observed in ID disordered soft-core boson system (^^=0.625) 
by varying the strength of on-site repulsion from zero (free) to strong (hard-core) couplings 

in- 

We made a semi-quantitative comparison of these numerical results with several results 
related to high-T^ superconductors. Hereafter we assume bosons to have a linear size equal 
to the coherence length on the Cu02 plane, a = C,ab — lOA. Recently Matsuda et al. 
pO| carried out transport measurements which provided unequivocal evidence of the KT 
transition in a one-unit-cell(12A)-thick YBa2Cu307_x> (YBCO) fllm with a maximum value 
of Tkt = 30K. The three-dimensional interlayer coupling, thus, is not a necessary condition 
for flnite-Tc superconductivity, although the interlayer coupling plays an additional role in 
elevating the |^ . By comparing the experimental and numerically-determined values of 
Tkt in the clean half-fllled system {2ml^ab'k-QT]^'i /h^ ~ 0.4), we estimate the boson mass 
as ml = 2m*j^j.j.igj, ~ 5.9me, which is comparable to the value of "^carrier — 2.5me given by 
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Schlesinger et al. from infrared measurements of YBCO. 

In high-Tc superconductors, chemical doping does not only result in carrier-density mod- 
ulation, but also introduces some degree of microscopic disorder in the superconductive 
Cu02 plane, although the latter effect has not been well elucidated experimentally. One 
remarkable feature commonly observed in experimental data is that the relationship be- 
tween the transition temperature Tc and the hole density 5 (per unit cell) follows a common 
bell-type curve: Tc{5) rises above doSsct — 0.05, reaches a maximum at 5Tc(max) — 0.18, and 
falls down to zero at 5end — 0.3. In addition, the Tc-(5 curve is symmetric about 5 = ^^^(maa:), 
possibly implying some hidden symmetry. These experimental features strongly correlate 
with the results presented in Fig. 2. Furthermore we theoretically obtained the universal 
parabolic curve, similar to the ansatz proposed in Ref. 0. By simply assuming that the 
boson particle density scales with hole density as rib = (^/2)(^afe/^af))^ with a unit cell size 
lab — 4A, we get Hb ~0.16, 0.56, and 0.94 for ^offset, ^Tcimax), and 6cnd, respectively. This is 
also consistent with the computational results. Thus, the existence of parabolic symmetry 
and the consistency of characteristic values of ml and 6 strongly suggest that superconduc- 
tivity within individual Cu02 layers is caused by 2D superfluidity of local bosons moving in 
somewhat disordered media. 

To summarize, our quantum Monte Carlo simulations have provided further support that 
the 2D hard-core bose gas undergoes a KT transition to become superfluid, and that the 
low temperature superfluid density shows a parabolic-like dependence on particle density. 
In disordered systems, we have shown that the bosons localize differently in the dense and 
dilute cases. The computational results correlate with the vs. 6 relationship for high-T^ 
superconductors and may partially illustrate the nature of the superconducting transition. 

We thank S. Nakajima, T. Ohtsuki, T. Uda, Y. A. Ono, and M. Hirao for their dis- 
cussions. T.O. thanks M. Suzuki, S. Miyashita, and M. Marcu for their comments from 
theoretical viewpoint, and Y. Matsuda for showing his experimental results prior to publi- 
cation. 
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FIGURES 

FIG. 1. (a) two-dimensional superfluid density Ug, (b) temperature derivative of Ug, and (c) 
specific heat Cv, as a function of temperature T, for a clean half- filled system. The symbols denote 
the cases for different lattice sizes L = 4 (circles), 6 (triangles), and 8 (squares). The solid curves 
represent spline fits to the calculated data. The dotted line in (a) corresponds to the universal 
jump condition of the Kosterlitz-Thouless transition. 

FIG. 2. Low-temperature superfluid density as a function of boson particle density, for different 
degrees of disorder. This simulation was performed at low temperature T = 0.25 t/kB on an 8 x 8 
lattice. Each value denotes the thermal and sample averages over five different disordered states. 
The degrees of disorder are denoted by circles for the clean system (A = 0), triangles (A = 1.0), 
squares (A = 1.25), and crosses (A = 1.75). The solid lines represent parabolic fits to the calculated 
data. 

FIG. 3. Localization in dense (n^ = 1/2) and dilute {ni, = 1/8) bose systems. The curves show 
the low temperature superfluid density as a function of disorder strength, at a fixed temperature 
T = 0.25 t/kB- The symbols correspond to lattice sizes L = 4 (black circles), 8 (black squares) for 
the dense system, and L = 4 (circles), 6 (triangles), 8 (squares) for the dilute system. The two 
insets represent the density distribution of bosons on the 8x8 lattice. The upper and lower insets 
correspond to the cases of (n^ = 1/8, A = 1.0) and (n^ = 1/2, A = 2.5), respectively. 
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